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Abstract 

We find exact solutions to the Dirac equation in D-dimensional de Sitter space- 
time. Using these solutions we analytically calculate the de Sitter quasinormal (QN) 
frequencies of the Dirac field. For the massive Dirac field this computation is similar 
to that previously published for massive fields of half-integer spin moving in four 
dimensions. However to calculate the QN frequencies of the massless Dirac field we 
must use distinct methods in odd and even dimensions, therefore the computation 
is different from that already known for other massless fields of integer spin. 
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1 Introduction 

The recent studies on the propagation of classical fields in D-dimensional (D > 4) space- 
times are mainly motivated by the Brane World scenario in String Theory and the scrutiny 
of the higher dimensional features of General Relativity. Also, the dS-CFT and AdS-CFT 
correspondences have motivated the study of classical fields moving in spacetimes with 
dimension different from four, generally in asymptotically de Sitter or anti-de Sitter back- 
grounds (see Ref. [1] for a review and for some references on this topic see [2]-[20]). 

*Electronic address: alfredo@fis.cinvestav.mx 
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As is well known, using classical fields as a probe we can explore the properties of a 
spacetime. An essential part of this analysis is the computation of the quasinormal (QN) 
frequencies, because they depend on the parameters that characterize the spacetime, for 
example, the mass, charge, and angular momentum [3]. The quasinormal modes (QNMs) 
were first computed for classical fields propagating in four dimensional black holes, because 
it is expected that the QN frequencies of the gravitational perturbations can be detected 
by the gravitational wave observatories and these are also useful to analyze the linear 
stability of the black holes [4]. 



Recently these modes have been helpful in other research lines as the AdS-CFT and 
dS-CFT correspondences of String theory and the thermodynamics of black holes [18] . 
[T9] . [20J. Also, as in four dimensions, the QNMs are useful tools to study the physical 
features of the D-dimensional backgrounds. These facts have motivated the computation 
of the QN frequencies for different fields moving in D-dimensional backgrounds [5J-[19|. 

A simple solution to the field equations of General Relativity is the de Sitter back- 
ground [2TJ. Owing to its simplicity the propagation of classical fields in this spacetime 
has been studied, see for example pj)] [IE], [21] [27J and references therein. Its QNMs 
are defined as the solutions to the equations of motion for the classical fields that satisfy 
the boundary conditions [17], [22], [23]: (i) the field is regular at r = 0; (ii) the field is 
purely outgoing near the cosmological horizon. 

The exactly solvable systems are usually limits of more realistic systems and allow us to 
study in detail some properties of a physical process and test some methods which can be 
used to analyze more complicated systems. Thus they are powerful tools in many research 
lines. Therefore we expect that the examples of exactly computed QN frequencies for D- 
dimensional spacetimes may play an important role in future research. For some examples 
of D-dimensional backgrounds whose QN frequencies have been exactly computed see 
Refs. P3HIE]. 

The exact computation of the QN frequencies for the de Sitter spacetime has attracted 
much attention and in Refs. [H]-[T8], [22] [25] have been exactly calculated for several 
fields and different spacetime dimensions. We expect that these values of the QN frequen- 
cies can be useful to explore the physical properties of the de Sitter spacetime. Our main 
objective here is to provide an additional field for which its de Sitter QN frequencies can 
be exactly calculated. 

Nevertheless, notice that in Refs. [H]-[T8], [22] [25] are reported different conclusions 
on the existence of these QN frequencies for the massless fields. For example, in Refs. 
[2"2"] . [24J is asserted that de Sitter QN frequencies of the massless scalar field do not 
exist in four dimensions. Also, Natario and Schiappa in Ref. [15] claim that for the 
gravitational perturbations the QN frequencies are well defined only in odd spacetime 
dimensions. In Refs. [T7J , [23] we showed that the de Sitter QN frequencies are well defined 
for electromagnetic and gravitational perturbations in D > 4 dimensions. Furthermore, 
in Ref. [25] Kim and Myung affirmed that the de Sitter QN frequencies do not exist; 
however in the latter reference they only search real frequencies and, as is well known, 
the QN frequencies are complex numbersQ 

For the massive scalar field in Ref. [16] Du et al. asserted that this field has well defined 

1 It is convenient to note that in Ref. [35] Kim and Myung impose slightly different boundary conditions 
on the fields. 
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de Sitter QN frequencies only if the mass of the field satisfy some conditions. Using a 
different definition for the radial flux of the scalar field to that of [16], we showed in Ref. 
[T7] that the QN frequencies are well defined without any restriction on the mass of the 
Klein-Gordon field whenever the test field approximation holds. Furthermore in Ref. [T8] 
were found well defined de Sitter QN frequencies for the scalar field in D-dimensions but 
Abdalla et al. did not analyze them in detail (see Ref. [17] for more details). 

As is well known, sometimes the Dirac field behaves in a different way than the integer 
spin fields when they propagate in a curved background. Thus for the Dirac field is 
interesting to compute its QN frequencies in de Sitter spacetime. In three and four 
dimensions these frequencies were already calculated in Refs. [33], [23]. Nevertheless, in 
the previous references for a Dirac field moving in D-dimensional de Sitter background, 
for D > 5, the de Sitter QN frequencies were not computed. 

Taking into account the recent results on the separation of variables for the Dirac 
equation in spherical symmetric Z)-dimensional backgrounds [28], [22], [30] we find exact 
solutions for this equation in D-dimensional de Sitter spacetime. Using these solutions 
we calculate the QN frequencies of the massive and massless Dirac field and therefore we 
extend the results of the previous references, filling a gap in the literature. We compute 
separately the QN frequencies for the massless and massive Dirac field in order to study 
in detail the former case, because for the massless fields there are different results on the 
existence of the de Sitter QN frequencies, as we already commented. 

We organize this paper as follows. In Sect. [2] we briefly describe the main results of Ref. 
[28] on the reduction of the Dirac equation to a system of ordinary differential equations 
when the D-dimensional background spacetime is spherically symmetric. In Sect. [3] we 
find exact solutions to the radial equations of the previous section in D-dimensional de 
Sitter spacetime and using them we calculate the QN frequencies of the Dirac field; we first 
study the massless case and then the massive case. In Sect. H]we make a brief discussion 
on the results obtained. 

In Appendix |A] we show that if the background spacetime is D-dimensional de Sitter 
spacetime, then the system of differential equations for the radial functions of the massless 
Dirac field given in Ref. [25] reduces to the system that we solve in Sect. [3J In Appendix 
IB1 we find the effective potentials of the Schrodinger type equations for the Dirac field 
moving in D-dimensional de Sitter spacetime and study some related facts. In Appendix 
ICl for the de Sitter QN frequencies of integer spin fields we explain and use an alternative 
method of computation. 

2 Dirac equation in D- dimensional de Sitter space- 
time 

Using a Cartesian gauge, Cotaescu in Ref. [28] proved that in D-dimensional spherically 
symmetric spacetime the Dirac equation reduces to a system of differential equations for 
the radial functions. In this section we make a brief summary of the results presented in 
that reference where more details can be consulted (see also [27]). 

To simplify the Dirac equation, in Ref. [28] Cotaescu write the metric of a D-dimen- 



3 



sional spherically symmetric spacetime in the formal 

ds 2 = W 2 dt 2 



w 2 vv 2 „ irv , 



w 



V 2 



where U, V, and W are functions of y, dQ 2 is the metric of the p-dimensional unit sphere 
and D = p + 2. Using the form for the metric of the background spacetime, in Ref. 
[28] is shown that the equation of motion for a Dirac field of mass M simplifies to 



MW 



-Uj- + k— 

ay y 



ay y 



-MW 



&{y) 




(2) 



if we use a Cartesian gauge and make the following ansatz for the Dirac field 

fi%)4> K {n p ) 



1 i, 
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(3) 



for even D, and 
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/^ M (2/)^NW+^l|(l/)0-|«|(^ 



(4) 



for odd D. In formulas ([3]) and (jlj), represents the angular spinors appropriate for 

the dimension of the spacetime under study and the quantity k is equal to [28J 



/, = ± ( ~ + n ) . 



n = 0,1,2, .... 



In static coordinates the D- dimensional de Sitter metric takes the form [21] 

- 2 ^ dr 2 



ds' 



r 



1-^ldt 2 



r 2 dtf p , 



(5) 



(6) 



where the quantity L is related to the cosmological constant. Thus the quantities U, V, 
and W for this metric are equal to [27], [28] 



V 



w 



1 



?/ Lsinh(|;) cosh(|;) 
Hence, taking into account the previous expressions, Eq. (T5]) simplifies to 
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2 In this manuscript we adopt a notation slightly different from that of Ref. 
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Employing the coordinate z = tanh(y/L), related to the static radial coordinate r by 
r = zL [21], we find that Eq. ([8]) transforms into 
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where uj = ujL and M 

From Eqs. ([9]) and defining the new functions R\ and R2 by 

we find that these functions satisfy 



(10) 
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This system of differential equations has the same mathematical form already found in 
Refs. [23], [26], where we studied the massive Dirac field moving in 3D and AD de Sitter 
spacetime. 

Another method for simplifying the massless Dirac equation to a system of differential 
equations for the radial functions is given in Ref. [29] (see also [SD])- The massive Dirac 
equation is not explicitly covered in Ref. [29]; therefore here we follow |28j. In Appendix IA1 
we show that if the background spacetime is the D- dimensional de Sitter, then the system 
of ordinary differential equations found by Das, et al. in Ref. [29] can be transformed into 
Eqs. dHJ> with M = 0. 



3 Quasinormal frequencies of the Dirac field 

In this section we find exact solutions to the system of differential equations ffTTT) and using 
them we calculate the QN frequencies of the Dirac field propagating in D-dimensional de 
Sitter spacetime. We first find these QN frequencies in the massless case and then in 
the massive case. For the massless and massive Dirac equations we use different methods 
to make the simplification of the radial equations to hypergeometric type differential 
equations. The procedure for the massless Dirac field is simpler than that for the massive 
Dirac field, and we must use a more complex method to find the exact solutions in the 
last case. 

3.1 Massless Dirac field 

Following a similar procedure to that used in Sect. 2 of Ref. [23] (see also [26]) we can find 
in a straightforward way the solutions to the system of differential equations (fTTj) when 
M is equal to zero. 
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As in Ref. [26], we make the ansatz for the radial functions R x and R 2 

R x =v Bl (l-v) Cl R x , R 2 = v B %l-vf 2 R 2 , (12) 



where 



2z 
l + z : 




iui ( lui 

' 2 ' 2 ' 



C 2 = I (13) 

iw _i_ 1 1_ iu)_ 

2 + 2' I 2 2 ' 

to find that the functions R x and R 2 must be solutions of the hypergeometric differential 
equations 

v(l - v)^f + (cj - (aj + b T + l)v)^- - ajbjRj = 0, (14) 

where / = 1,2, and the quantities eij, bj, and Cj are equal to 

„ 1 iu) „ 1 ids 

a x = B x + d + - - —, a 2 = B 2 + C 2 + - + — , 

b x = B x + C x + —, b 2 = B 2 + C 2 - —, 

ci = 2B X + 1, c 2 = 2B 2 + 1. (15) 

It is convenient to note that these solutions reduce to those already given in Refs. [23] , 
[26J for the massless Dirac field propagating in 3D and 4D de Sitter spacetime. 

As noted in Refs. [17], [23], [26] there are several possible values for the quantities a±, 
bx, and c x (a 2 , b 2 , and c 2 ) depending on the chosen values for B x and C\ [B 2 and C 2 ). 
In the following we only study the case B\ — | + n and C x — y + | (-B2 = f + n anc ^ 
C2 = -jf)- In the other cases we expect to find identical physical results. 

We first study the radial function R\. As in Refs. [T7J, [23], [26] we can show that the 
only solution of the differential equation ffT4j) . with 1 = 1, that leads to a radial function 
regular at r = is the hypergeometric function 

Ri = 2 F x (a x ,b x ;c x ;v). (16) 

The other solution to the hypergeometric differential equation leads to a radial function 
divergent at r = (see Refs. [T7], [23], [2S])- Thus 

uD 1 

i?i = ^ n+p/2 (l - v)"2" + 2 2 Fi(a!, 6 l5 c i; u), (17) 

is regular at r = 0. 

As is well known, when c x — a x — b x is not an integer the hypergeometric function 
satisfies 



t? ( u \ r(ci)r(d - ai - 6 X ) .A.i 1 

2 Fi(ai, 61; ci; u) = — — — 2-^1 (ai, h; a x + b x + 1 - c x ; 1 

r(ci - ai)r(ci - 61) 
r(c 1 )r(a 1 + & 1 - Cl ) _ , ci _ ai _ 6l 
r(a 1 )r(6 1 ) 1 ] 

x iF\{c\ — a x , c x - b x ; c x + 1 - a x - b x ; 1 — v), 
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where T(z) denotes the Gamma function. Exploiting this property of the hypergeometric 
function we can write the radial function (fITI) as 

[ iui 1 

R % = V n +P/ 2 Ul-v)2 +i 2 F 1 (a 1 , h; a x + h + 1 - c x ; 1 - v) 

r(ci)r(ci - ai - bi) r(ci)r(ai + &i - ci) 
x r( Cl - ai )r(ci - 60 r(oi)r(6i) 

iui 1 | 

x(l - 2 ~4 2 Fi(ci - ai,c a - h; cx + 1 - a x - bx, 1 - v) > . (19) 

In the previous expression the first term in curly brackets represents an ingoing wave 
near the cosmological horizon, while the second term represents an outgoing wave. Hence, 
to get a purely outgoing wave near the cosmological horizon we must satisfy the condition 

ci - di = — n 2 , or ci - bi = -n 2 , n 2 = 0,1,2, (20) 

We cannot fulfill the first condition in Eqs. ( 120|) . but from the second condition we obtain 
that the QN frequencies of the massless Dirac field moving in D-dimensional de Sitter 
background are equal tc(H 

p + 1 

iu = n-\ hn 2 . (21) 

Now, we note that for these frequencies 

c x - a x - h = - (n + 1 + | + n 2 ) , (22) 

hence c\ — ax — bx is not an integer for odd p, but for even p this quantity is an integer. 
So for even p the previous procedure does not produce the de Sitter QN frequencies of 
the massless Dirac field since in this case, if we impose the second condition of Eqs. (|20|) . 
then we contradict the assumption that the quantity c\ — ax — bx is not an integer already 
used to expand the radial function in formula (fT9j) [TT] . [23] . 

In Ref. [23] we noted that for a massless Dirac field moving in 3D de Sitter spacetime 
the QNMs boundary conditions can be satisfied when cx — ax — bx is not an integer. Thus 
our result agrees with that of the previous reference. 

For the massless Klein-Gordon field, the electromagnetic and gravitational perturba- 
tions we already showed that in even and odd dimensional spacetimes we cannot fulfill 
the boundary conditions for the de Sitter QNMs when the quantity cx — ax — bx is not an 
integer (see Sects. 2 and 3 of Ref. [17] and Sect. 3.2 and Appendix B of Ref. [23]). For the 
massless Dirac field we cannot satisfy the boundary conditions when c\ — ax — bx is not 
integer only in even dimensions. Therefore it shows a different behavior than the integer 
spin fields. 

As in Refs. [T7J, [23], for even dimensional de Sitter spacetimes we must analyze 
the case ci — ax — bx = —nx, nx = 1,2, 3,... Jj For these values of ax, bx, and cx the 

3 As in Refs. [17], [23] we also call to the quantities iu> the de Sitter QN frequencies. 
4 For other integer values of ci — a\ — b\ we do not get de Sitter QN frequencies; see Refs. [17], [23] 
for more details. 
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hypergeometric function satisfies [31] 



j-i/ j x r(ai + h - ni)r(m) . 

2 Fx(ai, 61 a x + bi -ni,v) = — (1 - v ) 

r(ai)r(6i) 



ni — 1 

X 

s=0 



E (ai - ni) s (6i - wi) s , (-l) wi r(ai + fei - ni) 

s!(l-n x ) s 1 vj rC^-nOr^-m) 
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00 

x 



s=0 



+ ^(ai + s) + tpih + s)], (23) 

where ip( z ) — d\nT(z)/dz, (z) = 1, and (z) s = z(z + 1) • • • (z + s — 1) for s > 1. From 
the previous formula and the radial function ffTTj) . to get a purely outgoing wave near the 
cosmological horizon we must impose the condition [T7|, [23] 

a\ — rix — — n 2 , or 61 — ni = — n 2 . (24) 

We cannot fulfill the second condition in Eqs. (1241) . but from the first condition we find 
that the QN frequencies of the massless Dirac field moving in D-dimensional de Sitter 
spacetime, for even D, are equal to 

p + 1 

iuj = n H h n 2 . (25) 



For the massless Dirac field propagating in AD de Sitter spacetime in Sect. 3 of Ref. 
we show that the QN frequencies are well defined. We also prove that the QN frequencies 
can be calculated when the quantity c\ — a\ — b\ is a negative integer. Hence our results 
agree with those for the massless Dirac field given in the previous reference. 

Thus from formulas (|2ip and (|25|) we can assert that for even and odd D the de Sitter 
QN frequencies of the massless Dirac field are well defined. Analyzing the radial function 
i?2, in even and odd dimensions, we also find the de Sitter QN frequencies (|25|) (or (I2~T|) ). 
It is convenient to notice that the quantities iu) of formula (|25|) are integers for odd p 
and half-integers for even p, as the de Sitter QN frequencies of a massless Klein-Gordon 
coupled to curvature with coupling constant £ = (D—l)/ {AD) (see expression (55) of [T7]). 
Nevertheless, if the mode number n 2 changes by the same quantity, then the increments 
in the imaginary part of the QN frequencies u calculated here for the massless Dirac field 
and in Ref. [17] for the massless scalar field are different. 

Notice that the QN frequencies fl25l) reduce to those calculated in Sects. 2 and 3 of 
Ref. [23J for the massless Dirac field moving in 3D and AD de Sitter spacetime (see also 
[TB]). We also point out that for two different odd dimensional de Sitter spacetimes the 
set of QN frequencies for the massless Dirac field are equal except for the first modes. 
This also happens for even dimensional de Sitter spacetimes. 



3.2 Massive Dirac field 

Exploiting a similar procedure to that given in Sect. 4 of Ref. [23j we can solve the system 
of differential equations (TTTT) when M^O. Thus in this section we find exact solutions to 
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the equations of motion for a massive Dirac field propagating in D-dimensional de Sitter 
spacetime and we use them to calculate its QN frequencies. 

Making the following ansatz for the radial functions Ri and R 2 [23J 



Ri = (l-2?)-i(l-zMf 1 -h), 
# 2 = (l-z 2 )-i(l + ^(/ 1 + / 2 ), 

and taking fi and / 2 in the form 



fi 
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Ox, 



where x = z 2 and 
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c ^l-x) B2 f 2 
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(27) 



'tUl I 1 

2 4' 
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(28) 



we find that the functions fi and / 2 are solutions of the hypergeometric differential equa- 
tions 



a;(l - x)—^ + (c/ - (a/ + 6/ + l)x)— - - a/6/// = 0, 
dar ax 



(29) 



where the quantities a/, 6/, and cj are equal to 

1 



ai = Si + d + 
h = B l + d + 



1 1 

4 + 2 

1 1 
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+ zM 



5 2 + C 2 + - + 



ci = 2Ci + 



1 



a 2 
62 

c 2 = 2C 2 + 



2 \2 



s 2 + c 2 + 
1 



1 
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-iM 



%M 



(30) 



Depending on the chosen values for the quantities Bi and C/, we get several values 
for aj, bf, and cj as already noted in Refs. [17], [23]. In the following we study the case 
B x — y + Ci — I , B 2 — y + j, and C 2 = ^ + f • 

The solutions of Eqs. fl29|) that lead to radial functions regular at r = are 

// = 2 F 1 (a I ,b I ;c I ;x). 



(31) 

The another linearly independent solution of Eq. (129]) leads to divergent radial functions 
at r = 0. Furthermore, from Eqs. ffTTl) we can show that if we choose fi as in formula 
d3U> then [23] 

A = — -2-F\(a 2 , b 2 ; c 2 ; x), (32) 

Cl 



and the radial function Ri is equal to 
R 1 = (l-z 



1 1 

"4(1 - z)-2 Z Z 



2 Fi(a 1 , bi, ci] z 2 ) H 2 Fi(a 2 , b 2 ; c 2 ; z 2 

Cl 



(33) 
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(a similar expression holds for the function R 2 ). 

Using the property of the hypergeometric function given in formula (ITBl) we can write 
the radial function Ri (1331 in a similar form to the complicated formula (92) of Ref. [23] . 
From this formula, to obtain a purely outgoing wave near the cosmological horizon we 
must satisfy the condition 

c 1 — a 1 + l — —n 2 , or c 1 — b 1 — —n 2 , n 2 = 0,1,2, (34) 

From the previous conditions and expressions fl30|) . we find that the QN frequencies for 
the massive Dirac field moving in Z)-dimensional de Sitter spacetime are equal to 

T) ~\~ 3 ~ T) ~\~ 1 ~ 

iu> = n + h 2n 2 - iM, iu = n + h 2n 2 + iM. (35) 

Applying a similar method we can analyze the radial function R 2 ; in this case we also 
find the de Sitter QN frequencies ( l35l) . Taking the limit M — > in expressions ( l35l) we get 
the frequencies for the massless Dirac field given in formula ( 1251) . It is convenient to note 
that for the QN frequencies f l25l) and ( 1351) the amplitude of the field diminishes with time. 
Also notice that the QN frequencies (1351) have the same mathematical form than those 
for a coupled to curvature massive scalar field with coupling constant £ = (D — 1)/(4D) 
which appear in formula (54) of Ref. |T7] . 

In AD de Sitter spacetime the QN frequencies (155]) take the form 

5 3 
iu = n + - + 2n 2 - iM, iuj = n + - + 2n 2 + iM. (36) 

The second set of quantities in formula (1361) is equal to the second set of QN frequencies 
given in expression (94) of Ref. [23J. To compare the other set of de Sitter QN frequencies, 
notice that in formula (94) of Ref. [23] we wrote the first set of the de Sitter QN frequencies 

as 

iu = j + 2m - iM, m = 0,1,2,..., (37) 

(j is a half-integer, j > |), which is wrong, because in the previous reference we had 
imposed the condition (first condition in Eqs. (93) of Ref. [23] ) 

c\-a,\ = -n 2 , (38) 

instead of the first condition in Eqs. (1341) . Hence the correct expression for these QN 
frequencies is 

iu = j + 2 + 2m - iM. (39) 

The QN frequencies ( 1391) coincide with those given in the first set of quantities in 
formulas (1361) . Furthermore, in Ref. [23] we asserted that in 4D de Sitter spacetime for 
each j, the frequencies iu = j — iM are QN. From the corrected expression (1391) for them, 
we infer that this affirmation is not true. 
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4 Discussion 



In this paper we found exact solutions to the radial equations of the Dirac field in D- 
dimensional de Sitter spacetime and using them we calculated its QN frequencies extend- 
ing the results of Refs. [16], |23j. For the massless case we get well defined expressions for 
the QN frequencies in even and odd spacetime dimensions. 

Moreover, our results for the de Sitter QN frequencies in the massless and massive cases 
reduce to those already calculated in Ref. [IT] for three and four dimensional spacetimes 
(see also Ref. [23]). Notice that for D > 5, the exact solutions to the Dirac equation and 
the values of the QN frequencies fl25l) and fl35|) have not been previously computed. 

Using two different sets of solutions to the Dirac equation for the massless and massive 
cases, we find the exact values of the de Sitter QN frequencies. For the massless Dirac field, 
we note that in odd dimensions the analysis is straightforward, while in even dimensional 
de Sitter spacetime, due to the non-standard arguments of the hypergeometric functions 
we must use more complex properties of these functions to calculate the de Sitter QNMs. 
Hence the massless Dirac field shows a different behavior than the integer spin fields 
already studied. 

Observe that taking the limit M — > of the result obtained for the QN frequencies of 
the massive Dirac field we find the expression corresponding to the massless field. Also, 
notice that the result for the massless field was calculated employing a different set of 
exact solutions. 

From the results of Refs. [17], [23], [26], and this paper we can deduce that the mass- 
less Klein-Gordon field, the massless Dirac field, the electromagnetic (EM) field and the 
gravitational perturbations have well defined QN frequencies when they propagate in D- 
dimensional de Sitter spacetime. We also notice that for these massless fields the de Sitter 
QN frequencies are purely imaginary. For massive fields they can be complex. For other 
spacetimes in which the QN frequencies of some fields are purely imaginary see Refs. [17] , 
[23], [32]. 

Summarizing the findings of Refs. [T7], [23], [26] and the previous sections we present 
the QN frequencies u for the D-dimensional de Sitter spacetime in Table 1. In this table 
n is an integer, n% — 0, 1, 2, . . . , M — ML, M being the mass of the field, q — 0, 1, and 2 
for tensor, vector, and scalar gravitational perturbations, respectively. 

To prove or disprove our results it would be interesting to obtain numerical values 
for the de Sitter QN frequencies. To our knowledge, for the Rarita-Schwinger and Proca 
fields propagating in D-dimensional de Sitter background, a similar analysis to that of 
Refs. [TTJ, [23], [26], and the present paper is missing. 

Is there any relation between the de Sitter QNMs and the thermodynamics properties 
of the cosmological horizon? We believe that this fact should be investigated additionally. 
For the case of a static black hole, see Refs. [2D] for more details. 
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Table 1: de Sitter QN frequencies 



r leld 


QN frequency u 






Massive scalar field 


-i L + 2n 2 + £±i + [(2±i) 2 


-M 2 ] l lA 


/L 


(n > 0) 


-< (n + 2n 2 + H±i-[(H±i) 2 


— M 2 ] l /j 


/L 


1\ T ' T \ • (2 11 

Massive Dirac field 
(n > 0) 


—i(n + ^ + 2n 2 — iM)/L 
-i(n + 2±i + 2n 2 + iM)/L 






T — 1 TV T C* 1 1 1 T 

EM field modes I 
(n > 1) 


-i(n + p- f + 2n 2 )/L 
-i(n + 2 + 2n 2 )/L 






EM field modes II 
(n > 1) 


—i(n + p + 2n 2 )/L 
-i[n + I + 2n 2 ) / L 






Gravitational perturbations 
(n > 2) 


—i(n +p + 1 - q + 2n 2 )/L 
—i(n + q + 2n 2 )/L 







A Das-Gibbons-Mathur equations 



The metric of a D-dimensional spherically symmetric spacetime in isotropic coordinates 
takes the form@ 

ds 2 = -f{f) dt 2 + g(f) [df 2 + f 2 dfij], (40) 
where the variable f is called isotropic radial coordinate and D = p + 2 as in Sects. [2] and 

For a D-dimensional spherically symmetric background whose metric is written using 
the isotropic coordinates of formula (HOI in Ref. |29] Das, et al. showed that the equations 
of motion for a massless Dirac field reduce to the system of ordinary differential equations 
(Eqs. (24) of Ref. 



h(f) 



d p + n 
df f 



h(r) 



d_ 

df 



Gi = iujF\ 



F} = iuGl 



(41) 

. For more details see 
(42) 

r 

that exists between isotropic and static coordinates for the de Sitter spacetime, we find 

•(-■••+) for the spacetime 



5 For consistency with Ref. [29] in this Appendix we use the signature ( 
metric. 



where h = \J f/g, F^, and G\ are functions of r and n = 0, 1, 2, 
Ref. |22|. 

Using the relation 

L + VL 2 - r 2 
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that the system of ordinary differential equations (j4Tj) becomes 



i_^)i_i_E±n 

dz z 
. 2 i d n 
dz z 



to-' 



;i -z 2 )^ 



-F 1 



-r 1 

i "n' 



where r = zL and a; = cjL as in Sect. [2j 

Replacing in Eqs. (j4"3"j) — > — G n , i 7 ^ — > F n , and taking 



G n — HG r 



where 



H 



i + VT 



-p/2 



we get that the functions G n and F n satisfy 



i d 



p\ 1 



2/ z 



1-^)5—- ( // -h - I - I (V 



;i -z 2 )t 

Defining the new radial functions 

Ri — G n + F n , R 2 = G n — F n} 

we obtain that the system of differential equations (I46p becomes 

id iuj \ „ f p\ 1 



-G n . 



dz i-z 2 U v 



dz 



'1-Z 2 )2 



i?9 



2/ z 
p\ 1 



i?2 



(43) 



(44) 
(45) 



(46) 



(47) 



Thus in D-dimensional de Sitter background, the system of differential equations 
the massless Dirac field simplify to Eqs. (1111) with M = and k — § + n. We 
solved this system of differential equations in Sect. [3j 



(48) 

(HI]) for 

already 



B Effective potentials 

For several spacetimes in which exact solutions to the equations of motion for the classical 
fields can be found, in Appendix C of Ref. [17] we noted that the effective potentials of 
the Schrodinger type equations are of the Poschl- Teller fornj§ 

^ ^ sinh 2 (x) cosh 2 (x) ^ 



6 Here the variable x is different from that of Sect. [31 
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We presented the values of A and B for these spacetimes in Table 1 of Ref. [Tj 

Applying the method described in Ref. [33], we can prove that Eqs. ( JTTj) for the 
massive Dirac field in de Sitter spacetime reduce to a pair of Schrodinger type equations 
with effective potentials 

14 = »' 2 ±^, (50) 
dr* 



(1- z 2 )^[k 2 + (Mzf}l 
z[k 2 + (Mz) 2 ] + ^lp!l' 

1 f M z s 

r* H arctan 



where 

W 



2u \ K 

«=p (51) 

z is defined in Sect. [21 r* is the usual tortoise coordinate for the D-dimensional de Sitter 
spacetime and the quantity W is known as a superpotential for the potentials V± 
For the massless field, we can write the effective potentials fl50|) as 



k 2 _ Kcosh(^) 
V± = sinh 2 (^) T Lsinh 2 (^)' (52) 

Therefore for a massless Dirac field moving in D-dimensional de Sitter spacetime, in these 
variables the effective potentials fl52l) are not of Poschl- Teller type (jUJ). Nevertheless, we 
note that the potentials fl52l) are of Rosen- Morse type (see Table I in Ref. [M]). This fact 
has not been previously observed^ 
Furthermore, the quantity 

W = k { ' =^csch(^), (53) 
z \L / 

is a superpotential for the potentials ( l52l) . Thus we can expect that the potentials ( !52l) 
have the same spectrum |34j . as we already noted. 

In Ref. p2] we showed that the modes I and II of the electromagnetic field have a 
different spectrum of the de Sitter QN frequencies (for fixed angular momentum number 
the QN frequencies have different parity in odd dimensions, while in even dimensions a 
finite number of QN frequencies, which depend on the spacetime dimension, are different). 



7 For the massless Dirac field propagating in the BTZ black hole we can simplify the equations of 
motion to Schrodinger type equations with effective potentials equal to (see expression (74) in [17) and 
Ref. [55]) 



V ± = m r* 



which arc of Morse type |34| . Furthermore, observe that the Poschl- Teller, Rosen-Morse, and Morse 
potentials are all shape invariant [34 . 
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For the electromagnetic field, the effective potentials of the modes I and II in D- 
dimensional de Sitter spacetime are equal to (Table 1 of Ref. [TT] ) 

1(1 4- n r>\ , (Jt^KD-A). (g-4)(J?-6) 

Vi = [ } 4 - 4 (54) 

sinh (x) cosh (x) ' 

Vn = l(l + D-3) + &=*^ _ (55) 
sinh 2 (x) cosh 2 (x) ' 

with x = r*/L. Hence we can write the Schrodinger type equation in the form [17] 

d 2 



dx 2 



+ (Co 2 + tt)- V 1 , ll $i,ii = 0, (56) 



where 

V x ,u = V I>n + Q, (57) 

and 

tt = (D -4 + Z) 2 . (58) 
Furthermore, we find that Vn can be written in the form (see formula (150p ) 

V-n = -^ + ^, (59) 



where 



n = ~ 2 J tanh(x) + ^ + Z - coth(x). (60) 



Therefore the SUSY related potential to Vh is 
Vu ~ dx +Wn 

(I 1)(Z I D 1) I ( J - 2 )p- 4 ) (g-4)(J-6) 

sinh 2 (a;) cosh 2 (x) ^ ' ^ 

which is not equal to V\. 

Thus the potentials Vi and Vn are not supersymmetric partner potentials. An inter- 
esting question is whether the modes I and II have a different spectrum of the de Sitter 
QN frequencies because their effective potentials are not supersymmetric partners. 



C de Sitter QNMs, an alternative method of compu- 
tation 

In Refs. [UJ, [23] we calculate the de Sitter QN frequencies of several massless fields 
with integer spin. Here we present a slightly different method of computation for these 
frequencies. We shall study in detail the case of the vector gravitational perturbations 
propagating in D-dimensional de Sitter background and at the end of this section we 
comment on the results obtained for the other massless fields. 
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In de Sitter spacetime the radial function R(z) for the vector type gravitational per- 
turbations satisfies the differential equation [15], [T7] 



where A = p(p — 2)/4, B = (2n + p — 2)/2, n an integer > 2, and z was already defined. 
Following Leaver [36] we take -R(z) in the form 

R{z) = z B+1 (l - z 2 )- & ' 2 R x {z) } (63) 

that fulfills the boundary conditions for the de Sitter QNMs if the function Ri (z) satisfies 
some conditions (see below), since z B+1 is a regular function at z = and (1 — z 2 )~ tu} / 2 
represents an outgoing wave near the cosmological horizon. Substituting expression ( l63i) 
into Eq. ( 162]) we find that the function Ri(z) must be a solution to 

d 2 R 1 (z) | / 1-ia) 2B + 2 dRi(z) + f 3iu + 2iuB + A + u 2 



dz 2 \ z + 1 z 2-1/ dz V 2(z + l 

2 + 3B + B 2 3B + B 2 - 3iu - 2iuB - A - u 2 + 2 



2(z + l) 2(z-l) 



i?i(»=0. (64) 



Making the change of variable x = 1 — z 2 we get that the function i?x (x) satisfies the 
ordinary differential equation 

,d 2 R 1 (x) / .„ / p 3\ \ diMx) 

x(l-x) — \l-iuj+\iuj-n---- 



dx 2 V. V 2 2/ / dx 



whose solution is [21] 

„ , v „ ,n-iu + p n+l-iu ._ 
i?i(s) =d 2 Fi( , ; 1 - ^; x) 

iCj ,n + iu +p n + 1 + iuj 
+ C 2 x luJ 2 F 1 ( , ;1 + iuj;x), (66) 

where C% and C 2 are constants. 

We note that the second term in formula (1661) produces an ingoing wave near the 
cosmological horizon, therefore we take C 2 = 0. Thus to satisfy the boundary conditions 
Ri(x) must be equal to 

R l (x) = C 12 F 1 ( , ;l-iu;x). (67) 

Also, notice that the parameters of the previous hypergeometric function satisfy 

c-a-6 = -- |-n. (68) 

Consequently the function Ri(x) is divergent at x = 1 (z — 0) and this divergence leads 
to a radial function divergent at the origin, except when R%(x) is a polynomial [31]. 
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Therefore in order that the radial function R(z) be regular at z = 0, we must satisfy the 
condition [3T] 



1 — iu + n 



n — iu + p 



712 = 0,1,2,... 



(69) 



2 



= -n 2 , 



or 



2 



= -7i 2 , 



Hence the de Sitter QN frequencies of the vector type gravitational perturbations take 
the form 



which are equal to those given in expression (37) of Ref. [TTj . Furthermore using the 
previous method we compute the de Sitter QN frequencies for the tensor and scalar 
gravitational perturbations, the electromagnetic modes I and II, and the massless Klein- 
Gordon field. The results that we get are the de Sitter QN frequencies already calculated in 
Ref. [IT]. For the massive Klein-Gordon field we also found the QN frequencies computed 
in the previous reference. 
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